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We study analytically the aging dynamics of the 0(n) model in the limit of 
n — > oo, with conserved and with non-conserved order parameter. While in the 
non-conserved dynamics, the autocorrelation function scales in the usual way 
C{t,tw) ~ C{t/tw), in the case of a conserved order parameter, 'multiscaling' 
manifests itself in the form C{t,tw) = C{h{t) / hit^)) , with a relaxation time 
growing more slowly than the age of the system (sub-aging), and h{t) a function 
growing faster than any length scale of the problem. In both cases, the effec- 
tive temperature associated to the violation of the fluctuation theorem tends to 
infinity in the asymptotic limit of large waiting times. 
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I. INTRODUCTION 

Aging of glassy systems is now well understood, at least from a qualitative point of view 
and different approaches have been used to understand such a behaviour. One of them 
is the interpretation of aging in terms of a coarsening process. The picture is the following: 
consider for instance an Ising ferromagnet, which is quenched at time t = Q below its critical 
temperature. When t increases, two types of domains emerge, with up and down spins. In 
the thermodynamic limit, equilibrium is never reached. At late times, domains have reached 
a typical size L{t). It is thus natural to assume scaling laws for the different quantities of 
interest For instance, one can try the ansatz S(k,t) ^ L'^g(kL) for the structure factor (in 
a d-dimensional space), or C{t,tw) ~ F{L(t) / L{tw)) for the two-time autocorrelation function, 
where g and F are scaling functions. The growth law L(t) determines then all the properties 
of the system. As an example, the droplet model for spin glasses |3j assumes a logarithmic 
growth, leading to C(t,t^) = F{\nt/ hitu,)- If the growth law is given by L{t) ~ P, like e.g. 
in a spinodal decomposition, one gets C{t,t') = F{t/tyj). This last behaviour is called 'simple 
aging' and has been analytically shown to hold within various non-random models 

Moreover, the above functional form for the correlation function is also found analytically in 
some mean-field models of spin glasses, which give the general form for the correlation functions 
in the aging regime C{t,tw) — C{-^^^), with h and C two scaling functions (valid in the 
two-time regime where both times are large, but with 1 < C < 0). Although the notations are 
different, the functional form is the same as in coarsening processes, and it is then very natural 
to try to interpret the ^.-function as a relevant length scale for spin glasses, as was done for 
instance in ref. Q . 

From the experimental and numerical side, it is found that a simple aging behaviour describes 
the data well, in many different systems. This is interpreted by saying that the relaxation time 
tr{tw) of the system scales as the age of the sample: tr ^ t^. However, a more subtle effect 
may appear, since very often grows more slowly than tw This effect has been called sub-aging 
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. In his pioneering experiments on polymer glasses, Struik |^] introduced the exponent ^ from 
the relation tr ~ t(i, with /i < 1. Different values of /i have been reported: Struik used /i ~ 0.9, 
experiments in spin glasses /i ^ 0.97 simulations of a structural glass were fitted using the 
value /X ^ 0.88 1^, and recently, experiments on a gel gave ^ ~ 0.9 [||. It can be checked (this 
point is discussed in detail in ref. [pj) that the /i-exponent is equivalent to the following choice 
of the /i-function: h{t) = exp(Y^(^)^~^). In accordance to what has been said above, this 
equivalence holds when ty^ ~* oo and t — ^ ^w- Another function, the 'enhanced power law' 
form h{t) = exp(ln"(t/to)) with a > 1, has been phenomenologically introduced in the context 
of spin glasses B, and the value a = 2.2 was used to fit experiments. This in turn gives the 
relation t^. ~ twT^'^'^~^ {^w) , valid in the regime — > oo and t — ^ / ln'^^^(t„,). 

These choices are nonetheless not clearly motivated from a theoretical point of view, since 
the mean- field spin glass models discussed above only predict the existence of h{t), and its 
analytical computation remains at present an open problem. In this context, simple models 
where h can be computed are much needed, but there are only few examples where sub-aging 
appears. Very recently, a model exhibiting sub-aging has been proposed by Rinn et al [|l0| , 
who studied a slight variation of Bouchaud's trap model for aging. This has given a theoretical 
support to the use of an exponent /i, even if its physical origin remains somewhat unclear. 
A scaling approach to the diffusion of a point particle in a low dimensional space has been 
proposed in ref. |pl]| , and leads in some cases to a sub-aging which can be well described by an 
enhanced power law. 

We study in the present paper a model for coarsening (the 0{n) model in the large-n limit) 
which also exhibits a sub-aging scaling in the autocorrelation function when the order parameter 
is not conserved. Its origin is the simultaneous presence in the system of two different length 
scales, whose consequence is the breakdown of the simple scaling laws generally used in domain 
growth processes. In particular, no t/<„,-scaling is found, and the relaxation time grows as tr ~ 
i„,/Vln (sub-aging). The autocorrelation is shown to be well represented in the asymptotic 
regime by an enhanced power law with a = 3/2, i.e. h{t) = exp((ln a;)"^/^). Interestingly enough, 
h{t) can not be interpreted in our example as a length scale. We do not want to argue that the 
model is a realistic one for the aging of polymers or spin glasses, but rather to give a possible 
physical explanation (the role of length scales [^) for the absence of the 'naive' i/i^-scaling, 
and exhibit a simple example where the /i-function can be computed and discussed in terms of 
length scales, which has not been done so far. 



II. THE 0{N) MODEL 

This model is one of the few exactly solvable model for coarsening. It was first studied by 
Coniglio and Zannetti [|l3| , who computed the scaling properties of the structure factor during 
the domain growth process. They pointed out the presence of the two mentioned length scales, 
and named 'multiscaling' the breakdown of the usual S{k,t) ~ L'^gikL). Bray and Humayun 
have shown, however, that this multiscaling was a peculiarity of the large-n limit, and proved 
that for a large but finite value of n, a 'normal scaling' was recovered ||l^. On an other hand, this 
'pathology' has been shown to appear as a relevant preasymptotic effect in different coarsening 
models like for instance the kinetic Ising model. 

The model is defined through the Hamiltonian 

H[cj>\ = I d'x (i(V0)2 + i-(n - 0^)^) , (2.1) 
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where (t>{x, t) is a n-component vector field in a d-dimensional space. Two different dynamics 
may be associated to this model, depending on whether or not the order parameter is conserved. 
In the case of a non-conserved order parameter, the dynamics is given by the so-called time 
dependent Ginzburg-Landau equation 



d<f){x,t) _ SH 



dt S4>(x,t) 



Vix,t), (2.2) 



where 'q{x,t) is a random Gaussian variable with mean zero and variance given by 
{'n{x, t)T]ix', t')) = 2T5{t - t')5{x - x'). 
For conserved fields, we add — in front of the r.h.s to get the Cahn-Hilliard equation 

^'^("'*)-V^('-^Ur,(.,t), (2.3) 



dt \S4>(x,t) 

where the variance of ri{x,t) is {T]{x,t)T]{x' ,t')) = —2TS{t — t')W^S{x — x'). 

We shall see below that the limit n ^ oo allows to solve the dynamics in both cases. The 
key point that makes the model exactly soluble is that in the limit of n ^ oo, the replacement 
(p^ /n (0^), where is one of the components of </>, can be made. The two types of dynamics 
are now successively considered. 



III. NON-CONSERVED ORDER PARAMETER: SIMPLE AGING 



The time dependent Ginzburg-Landau equation (2.2) associated to the Hamiltonian (2.1) is 



^ = VV + 0--(0')0 + »7, (3.1) 
ot n 

where the dependence on space and time has been removed for clarity. This differential equation 
is associated with random initial conditions, in order to reproduce the quench experiment 
described in the introduction, and (/)(a;, 0) is taken from a Gaussian distribution with zero 
mean and variance (</i(cc, 0)(/i(a;', 0)) = A^(a; — x'). From now on, we work at T = 0. In 
the coarsening problem, temperature does not play an essential role, provided it is below the 
critical temperature. The scaling regime can then be directly studied at T = 0. The review 
paper |^ provides a longer discussion of that point, and we discuss below how our results may 
be (slightly) changed by an non-zero temperature. 

The large-n limit results in the following equations which have to be self-consistently solved: 

|^=V20 + a(O0; a(t)-l-(02). (3.2) 

The solution is discussed in Refs. |p|,p^, and one finds for the Fourier transform (/)(fe,i) = 
/ dfix<f)[x,t)e-'''-=' 

(/.(fc,i) = 0(fc,O)e-'= M -j , (3.3) 

where io = A'^^^/Sir. It is now easy to compute the structure factor 

S{k,t) ^ ^{<t>{Km-k,t)) ^ {^i.ty/^e-^'^'K (3.4) 
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We used ((/)(fc, O)0(— fc, 0)) = AV^ from initial conditions. The structure factor may be written 
as S{k,t) = L'^g{kL), with L{t) = t^^^ and g{x) = {Sny/^ exp{-2x^), demonstrating the 
vahdity of the scahng hypothesis in that case. 
The autocorrelation function is defined as 



V 



d'^x{(f>{x,t)(f>{x,t^)) 



V J {2t:Y 



and may be easily computed: 



t + tu 



d/2 



Defining the scaling variable Ai = t/tw, C{t, t^) can be rewritten 



= i^i(Ai); Fi{x) 



1 + x 



(3.5) 



(3.6) 



(3.7) 



This last equation means that the autocorrelation function exhibits a simple aging behaviour. 
We have then illustrated on a concrete model the scaling approach to domain growth described 
in the introduction. We shall see in the next section the differences arising when sub-aging is 
present. 

Let us note here that a finite temperature does not affect the above discussion, since it 
simply introduces a short-time relaxation in the correlation function, that does not depend on 
the waiting time and corresponds to an equilibrium relaxation inside the growing domains. 
The long-time relaxation we are interested in, and which corresponds to the growth of the 
domains themselves is still described by (tJ.Tl). 



IV. CONSERVED ORDER PARAMETER: SUB-AGING 



The Cahn-Hilliard equation (2.3) associated to the Hamiltonian (2.f) is given by (still at 
T = 0) 



dt 



and is solved following the same steps as previously, leading to P,p^ 
m t) = m 0) exp (^-kh + M^)^ 
with to = A"'/''/(167r)^. The structure factor reads in that case 

s{k,t) ^ [L,{trY^''^^'^\ 



(4.1) 



(4.2) 



(4.3) 



where f{x) = 2x^ — x^. In this expression, two characteristic length scales have been defined: 
Li{t) = t^^^, and L^it) = (^-j^^jt^f ^ ^ ■ In the standard scaling form, S{k,t) ~ L'^g{kL), the 
structure factor varies as L'' with a prefactor depending on the scaling variable kL, whereas for 
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the multiscaling form ( |4.3| ), S varies as L", with an exponent a which depends continuously on 
the scahng variable kL2- The two scalings are thus completely different. 

Coniglio and Zannetti have interpreted this multiscaling in terms of domains composed 
of sub-domains, each sub-domain growing at a different rate. The initial motivation for the 
present work was indeed to investigate the possible existence of a 'hierarchy' of time scales, 
similar to the one found in mean- field spin glass models ('ultrametricity in time') [|l|,pi)[[l6[ . A 



different effect arises instead. Using eq.(4.2), one easily gets for the autocorrelation function 
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It is obvious from this expression that C{t^tw) cannot be written as a function of t/t^ only. 
The physical key ingredient for the absence of the usual scaling is the presence of two different 
length scales in the system. 



We prove now analytically that eq.(4.4) implies sub-aging. It has to be remarked first that 



when the time difference t = t — is equal to , one has 

In the asymptotic limit of large waiting times, the relaxation of C{t,tw) is complete in times 
T ^ tyj. In that regime, one can show that 

C{t,t^) ^^^^e^p\-^^^[-] ). (4.6) 




Defining the scaling variable A2 — r\/ln t^/t^, eq.(4.4) can finally be rewritten 



C{t,t^) ^ F,{X2); F2(a;)^exp(-— ). (4.7) 



The relaxation time grows hence as tr ^ t^/y/hi tw, i.e. more slowly than tw- this is a sub-aging 
behaviour. It is moreover possible to compute the function h{t) discussed in the introduction. 

The scaling form C{t,tw) — C (^^^^'^ should be valid in the two-time regime where both times 

are large, but with a non-zero value of the correlation function. In the present case, this regime 
is characterized by 

t^^oo, T~-^. (4.8) 

We have seen that a natural choice for h{t) would be Li{t) or L2{t), i.e. a length scale, since it 
is a common interpretation. This does not work, and a more complicated form has to be found. 
It is straightforward to realize that a possible choice is an enhanced power law: 

C{x)=exp(^~^\n'{x)y, ^t) = exp ((Int)^/^) . (4.9) 

The function h is neither Li nor L2, but a combination of the two, and therefore does not have 
a direct physical interpretation: h{t) ^ exp ( (Li/L2)^ 
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V. RESPONSE FUNCTIONS: INFINITE EFFECTIVE TEMPERATURES 



It is also relevant to study the response functions for aging systems, since it is a major 
prediction of the dynamical mean-field theory for spin glasses that interesting informations 
are encoded in the susceptibilities Up to now, we have studied aging in the two-time 

correlation functions C{t,tw)- In glassy systems, aging is also found in the related response 
functions R{t,tw), associated with a breakdown of the fluctuation dissipation theorem which 
at equilibrium would be TR(t,tw) = dC{t,tw) / dtw This is taken into account by introducing 
an effective temperature Toff through p7[ 



dC(t,t^) 

Tcsiq) = lim , 

tui^OC R[t^ tw) 



(5.1) 

C(t,t„)=g 



In coarsening systems, however, response functions have been shown numerically and ana- 
lytically to be weak, in the sense that Tgg oo . This property has been related to the 
decreasing density of topological defects (domain walls) during the coarsening. In the case of 
the 0{n) model, no topological defects are present if n > d, which is naturally the case in the 
large-n limit. We compute then R{t,tw) in the both cases studied above to obtain T^s. We 
refer the reader to ref. for the method, since we follow exactly the same steps. We get the 
two following expressions: 



R{t,t^)-^[-] [j^] . (5-2) 



in the non-conserved case, and 



(^/t \nt — -^tyj In tj\ 




in the conserved case (we dropped out all numerical constants). Combining eqs. (3. 6, 4. 4, 5.2,5.3), 
it is easy to show that for the non-conserved and the conserved case successive' 



y, one has: 



(<i-2)/4 



Teff(g) ~ lim t'^^-\ Teff(g) ^ lim Z --=-. (5.4) 

t™-+oo tu,^oo (lniu,)'.''+^''/* exp(vlnt„,) 

This holds for < g < 1, and shows that for d > 2, although there is no interpretation here 
in terms of defects, the effective temperature is infinite, as has been found so far in all domain 
growth processes |4| |l^ , [l9| . 

We studied in this letter the aging dynamics of the 0{n) model in the large-n limit. We 
showed that when the order parameter is not conserved, standard scaling laws hold, leading to 
a simple aging behaviour. We investigated the more interesting case of a conserved dynamics, 
and were able to show that the multiscaling observed in the structure factor does not imply 
a hierarchy of time scales ('ultrametricity in time' fl^). Rather, the relaxation takes place 
in a time scale which is shorter than the waiting time, tr ~ tu,/ln°~^(t^) with a = 3/2, the 
correlation function being well represented in that regime by C(t,tyj) = C {h(t) / hlt^)) , where 
h is an enhanced power law h{t) — exp(ln°(t)). This simple example exhibits then a very rich 
aging behaviour, whose origin is the presence of two different length scales during the coarsening 
process. It shows also that the interpretation of h{t) as a length scale may in some cases be 
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misleading. The enhanced exponential form that has been successfully used to fit spin glass 
experiments arises naturally from our computation. It implies that the relaxation time scales 
as tr ~ tyj / \n°'~^ (tyj) , which could hardly be experimentally distinguishable from a power law 
tr ^ tli^, when is very near to one, as it is in spin glasses. 
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